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Abstract. Gravitational waves detected from well-localized inspiraling binaries would al- 
low to determine, directly and independently, both binary luminosity and redshift. In this 
case, such systems could behave as "standard candles" providing an excellent probe of 
cosmic distances up to z < 0. 1 and thus complementing other indicators of cosmological 
distance ladder. 
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1. Introduction 

Coalescing binaries systems are usually con- 
sidered strong emitter of gravitational waves 
(GW), ripples of space-time due to the pres- 
ence of accelerated masses in analogy with 
the electromagnetic waves, due to acceler- 
ated charged. The coalescence of astrophysi- 
cal systems containing relativistic objects as 
neutron stars (NS), white dwarves (WD) and 
black holes (BH) constitute very standard 
GW sources which could be extremely use- 
ful for cosmological distance ladder if phys- 
ical features of GW emission are well deter- 
mined. These binaries systems, have compo- 
nents that are gradually inspiralling one over 
the other as the result of energy and an- 
gular momentum loss due to (also) gravita- 
tional radiation. As a consequence the GW fre- 



quency is increasing and, if observed, could 
constitute a "signature" for the whole sys- 
tem dynamics. The coalescence of a com- 
pact binary system is usually classified in 
three stages, which are not very well delim- 
ited one from another, namely the inspiral 
phase, the merger phase and the ring-down 
phase. Temporal interval between the inspiral 
phase and the merger one is called coalescing 
time, interesting for detectors as the American 
LIGO (Laser Interferometer Gravitational- 
Wave Observatory) (Abramovici et al. 1992) 
and French/Italian VIRGO (Caron et al. 1997). 
A remarkable fact about binary coalescence 
is that it can provide an absolute mea- 
surement of the source distance: this is 
an extremely important event in Astronomy 
(Capozziello et al. 2010). In fact, here we want 
to show that such systems could be used as re- 
liable standard candles. 
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2. Coalescing binaries as standard 
candles 

The fact that the binary coalescence can pro- 
vide an absolute measurement of the distance 
to the source, can be understood looking at 
the waveform of an inspiraling binary; as long 
as the system is not at cosmological distances 
(so that we can neglect the expansion of the 
Universe during the propagation of the wave 
from the source to the observer) the wave- 
form of the GW, to lowest order in v/c is (see 
(Maggiore 2007) for a detailed exposition of 
GW theory) 

ft+W= ^I(5^) 2/3 (i±-£i) cos[a , fe)] (1) 

ft* it) = fli (^) 2/3 cos i sin [<D (fe)] (2) 

where 3{ = 4(^) 5/3 , h + and h x are 
the amplitudes for the two polarizations of the 
GW, and i is the inclination of the orbit with 
respect to the line of sight, 



M, 



(mim 2 ) 3/5 



(3) 



(mi + m2)'/ 5 

is a combination of the masses of the two stars 
known as the chirp mass, and r is the distance 
to the source; / is the frequency of the GW, 
which evolves in time according to 



(4) 



t is retarded time, and the phase <t> is given by 



<D(f) = 2tt 



fdt'f(t'). 

Jin 



(5) 



For a binary at a cosmological distance, i.e. at 
redshift z, taking into account the propagation 
in a Friedmann-Robertson- Walker Universe, 
these equations are modified in a very simple 
way: 

1. The frequency that appears in the above 
formulae is the frequency measured by the 
observer, / b S , which is red-shifted with re- 
spect to the source frequency f s , i.e. / t, s = 
f s l{\ + z), and similarly t and t ret are mea- 
sured with the observer's clock. 



2. The chirp mass M c must be replaced by 
M c = (l+z)M c . 

3. The distance r to the source must be re- 
placed by the luminosity distance diiz). 

Then, the signal received by the observed from 
a binary inspiral at redshift z, when expressed 
in terms of the observer time t , is given by 

1 + cos^ i 

h + (?) = h c (t R ) cos [<D (t R )] , (6) 



h x (f) = h c (t R ) cos i sin [O (t R )] 



(7) 



where 

u (+\ 4 / GMg (z) n5 
he (0 = jTT 2 — 

d L (z) \ c 2 

Let us recall that the luminosity distance c/ L 
of a source is defined by 



(8) 



And 2 , 



(9) 



where T is the flux (energy per unit time per 
unit area) measured by the observer, and X 
is the absolute luminosity of the source, i.e. 
the power that it radiates in its rest frame. 
For small redshifts, di is related to the present 
value of the Hubble parameter Ho and to the 
deceleration parameter qo by 

^ =z+ l ( l_ rf + .... (10) 

c 2 

The first term of this expansion give just the 
Hubble law z - (Ho/c)di, which states that 
redshifts are proportional to distances. The 
term 0(z 2 ) is the correction to the linear law 
for moderate redshifts. For large redshifts, the 
Taylor series is no longer appropriate, and the 
whole expansion history of the Universe is en- 
coded in a function diiz). As an example, for a 
spatially flat Universe, one finds 



d L (z) = c(l+z) 



f 

Jo 



dz' 



(11) 



where H(z) is the value of the Hubble param- 
eter at redshift z. Knowing di(z) we can there- 
fore obtain H(z). This shows that the luminos- 
ity distance function d^z) is an extremely im- 
portant quantity, which encodes the whole ex- 
pansion history of the Universe. 
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Now we can understand why coalescing 
binaries are standard candles. Suppose that we 
can measure the amplitudes of both polariza- 
tions h+, h x , as well as / Q b s (for ground-based 
interferometers, this actually requires cor- 
relations between different detectors). The 
amplitude of h+ is h c (l + cos 2 i)/2, while the 
amplitude of h x is /z c cost. From their ratio 
we can therefore obtain the value of cos t, that 
is, the inclination of the orbit with respect to 
the line of sight. On the other hand, (4) (with 
the replacement M c — > M c mentioned abobe) 
shows that if we measure the value of / b s 
corresponding to a given value of f b s , we 
get M c - Now in the expression for h+ and h x 
all parameters have been fixed, except cIl{z). 1 
This means that, from the measured value 
of h+ (or of h x ) we can now read di- If, at 
the same time, we can measure the redshift z 
of the source, we have found a gravitational 
standard candle, and we can use it to measure 
the Hubble constant and, more generally, 
the evolution of the Universe (Schutz 1986). 
The difference between gravitational standard 
candles and the "traditional" standard candles 
is that the luminosity distance is directly 
linked to the GW polarization and there is no 
theoretical uncertainty on its determination a 
part the redshift evaluation. Several possibil- 
ities have been proposed. Among these there 
is the possibility to see an optical counterpart. 
In fact, it can be shown that observations 
of the GWs emitted by inspiralling binary 
compact systems can be a powerful probe 
at cosmological scales. In particular, short 
GRBs appear related to such systems and quite 
promising as potential GW standard sirens 
(Capozziello et al. 2010)(Dalal et al. 2006)). 
On the other hand, the redshift of the binary 
system can be associated to the barycenter of 
the host galaxy or the galaxy cluster as we are 
going to do here. 



1 It is important that the ellipticity of the orbit 
does not enter; it can in fact be shown that, by the 
time that the stars approach the coalescence stage, 
angular momentum losses have circularized the or- 
bit to great accuracy. 



3. Numerical results 

We have simulated several coalescing binary 
systems at redshifts z < 0.1 In this analysis, 
we do not consider systematic errors and errors 
on redshifts to obviate the absence of a com- 
plete catalogue of such systems. The choice of 
low redshifts is due to the observational limits 
of ground-based interferometers like VIRGO 
or LIGO. Some improvements are achieved, 
if we take into account the future generation 
of these interferometers as Advanced VIRGO 2 
and Advanced LIGO 3 . Advanced VIRGO is a 
major upgrade, with the goal of increasing the 
sensitivity by about one order of magnitude 
with respect to VIRGO in the whole detection 
band. Such a detector, with Advanced LIGO, is 
expected to see many events every year (from 
10s to 100s events/year). In the simulation pre- 
sented here, sources are slightly out of LIGO- 
VIRGO band but observable, in principle, with 
future interferometers. 

Here, we have used the redshifts taken by 
NED 4 (Abell 1989), and we have fixed the 
redshift using z at the barycenter of the host 
galaxy/cluster, and the binary chirp mass Mc, 
typically measured, from the Newtonian part 
of the signal at upward frequency sweep, to ~ 
0.04% for a NS/NS binary (Cutler et al. 1994; 
Capozziello et al. 2010). The distance to the 
binary di ("luminosity distance" at cosmolog- 
ical distances) can be inferred, from the ob- 
served waveforms, to a precision ~ 3/p < 
30%, where p = S /N is the amplitude signal- 
to-noise ratio in the total LIGO network (which 
must exceed about 8 in order that the false 
alarm rate be less than the threshold for de- 
tection). In this way, we have fixed the charac- 
teristic amplitude of GWs, and frequencies are 
tuned in a range compatible with such a fixed 
amplitude, then the error on distance luminos- 
ity is calculated by the error on the chirp mass 
with standard error propagation. 

The systems considered are NS-NS and 
BH-BH. For each of them, a particular fre- 
quency range and a characteristic amplitude 
(beside the chirp mass) are fixed. We start 

2 (?) 

3 (?) 

4 NASA/IPAC EXTRAGALACTIC DATABASE 
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with the analysis of NS-NS systems (Mq — 
\22M & ) with characteristic amplitude fixed 
to the value 1CT 22 . In Table 1, we report the 
redshift, the value of he and the frequency 
range of systems analyzed. In Fig. 1, the de- 
rived Hubble relation is reported. 



Object 


z 


K 


Freq. 

(Hz) 


NGC5128 


0.0011 


\<r' a 


0-=- 10 


NGC 1023 Gr. 


0.0015 


io- 22 


0-=- 10 


NGC 2997 


0.0018 


io- 22 


5 h- 15 


NGC 5457 


0.0019 


io- 22 


10 -=-20 


NGC 5033 


0.0037 


io- 22 


25 -=-35 


Virgo CI. 


0.0042 


io- 22 


30 -=-40 


Fornax CI. 


0.0044 


io- 22 


35 -=-45 


NGC 7582 


0.0050 


io- 22 


45 -=-55 


Ursa Major Gr. 


0.0057 


io- 22 


50 -=-60 


Eridanus CI. 


0.0066 


io- 22 


55 -=-65 



Table 1. Redshifts, characteristic amplitudes, 
frequency range for NS-NS systems. 




1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 



Fig. 1. Luminosity distance vs redshift for sim- 
ulated NS-NS systems. 

The Hubble constant value is 72 + 
1 km/sMpc in agreement with the recent 
WMAP estimation ( Larson, D., et.al. 201 1). 
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Object 


z 


h c 


Freq. 

(Hz) 


Pavo-Indus 


0.015 


io- J1 


65 -=-70 


Abell 569 


0.019 


io- 21 


75 -=- 80 


Coma 


0.023 


io- 21 


100 -=- 105 


Abell 634 


0.025 


io- 21 


110-=- 115 


Ophiuchus 


0.028 


io- 21 


130-=- 135 


Columba 


0.034 


io- 21 


200 -=- 205 


Hercules 


0.037 


io- 21 


205 -=-210 


Sculptor 


0.054 


io- 21 


340 -=- 345 


Pisces-Cetus 


0.063 


io- 21 


420 -=- 425 


Horologium 


0.067 


io- 21 


450 -=- 455 



Table 2. For each cluster we indicate redshifts, 
characteristic amplitudes, frequency range for 
BH-BH systems. 




0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 0.055 0.06 
Redshift 



Fig. 2. Luminosity distance vs redshift for sim- 
ulated BH-BH systems. 



The same procedure is adopted for BH-BH 
systems (M c = 8.67M , h c = IO -21 ). In 
Tables 2 we report the redshift, the value of 
he and the frequency range. The simulations 
are reported in Fig. 2, and the Hubble con- 
stant value computed by these systems is 69 + 2 
km I sMpc. 
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4. Conclusions 

We have considered simulated binary systems 
whose redshifts can be estimated consider- 
ing the barycenter of the host astrophysical 
system as galaxy, group of galaxies or clus- 
ter of galaxies. In such a way, the standard 
methods adopted to evaluate the cosmic dis- 
tances (e.g. Tully-Fisher or Faber-Jackson re- 
lations) can be considered as "priors" to fit the 
Hubble relation. We have simulated, for ex- 
ample, NS-NS, and BH-BH binary systems. 
Clearly, the leading parameter is the chirp mass 
M c , or its red-shifted counter-part M c , which 
is directly related to the GW amplitude. The 
adopted redshifts are in a well-tested range 
of scales and the Hubble constant value is 
in good agreement with WMAP estimation. 
The Hubble-luminosity-distance diagrams of 
the above simulations show the possibility to 
use the coalescing binary systems as distance 
indicators and, possibly, as standard candles. 
The limits of the method are, essentially, the 
measure of GW polarizations and redshifts. 
Besides, in order to improve the approach, 
a suitable catalogue of observed coalescing 
binary-systems is needed. This is the main dif- 
ficulty of the method since, being the coales- 
cence a transient phenomenon, it is very hard 
to detect and analyze the luminosity curves of 
these systems. Furthermore, a few simulated 
sources are out of the LIGO- VIRGO band. 

Next generation of interferometer (as 
LISA 5 or Advanced- VIRGO and LIGO) could 
play a decisive role to detect GWs from these 
systems. At the advanced level, one expects to 
detect at least tens NS-NS coalescing events 
per year, up to distances of order 2Gpc, mea- 
suring the chirp mass with a precision better 
than 0.1%. The masses of NSs are typically 
of order 1.4M Q . The most important issue that 
can be addressed with a measure of cIl(z) is 
to understand "dark energy", the quite myste- 
rious component of the energy budget of the 
Universe that manifests itself through an ac- 
celeration of the expansion of the Universe at 
high redshift. This has been observed, at z < 
1 .7, using Type la supernovae as standard can- 
dles (Riess et al. 1998; Perlmutter et al. 1999). 

3 (http://www.lisa-science.org) 



A possible concern in these determinations is 
the absence of a solid theoretical understanding 
of the source. After all, supernovae are com- 
plicated phenomena. In particular, one can be 
concerned about the possibility of an evolu- 
tion of the supernovae brightness with redshift, 
and of interstellar extinction in the host galaxy 
leading to unknown systematics. GW stan- 
dard candles could lead to completely inde- 
pendent determinations, and complement and 
increase the confidence of other standard can- 
dles, (Holz et al. 2005), as well as extending 
the result to higher redshifts. In the future, the 
problem of the redshift could be obviate find- 
ing an electromagnetic counterpart to the coa- 
lescence and short GRBs could play this role. 

In summary, this new type of cosmic dis- 
tance indicators could be considered comple- 
mentary to the traditional standard candles 
opening the doors to a self-consistent gravita- 
tional astronomy. 
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